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Abstract 

In  order  to  find  out  if  silicon  nanostructured  materials  resist  lithium  cations  insertion  and 
extraction  under  the  presence  of  germanium  atoms,  we  have  explored  theoretically  the 
potential  energy  surface  of  Sif,Gc\  Si6GeLi5\  Sir,GcLi7',  ShGcLi  |0\  SieGe^’,  Si6Ge2Li5',  and 
Si loGcioLi |(f  clusters.  This  study  was  performed  using  the  Kick  Coalescence  method 
coupled  to  Gaussian  09  program,  and  re-optimization  all  structures  with  the  B3LYP/def2- 
TZVPP  method.  Our  results  confirm  that  incorporating  germanium  atoms  to  silicon  clusters 
improves  and  prevents  fragmentation. 


Introduction 

Rechargeable  Li-ion  batteries  are  the  most  promising  energy  storage  devices  for  hybrid, 
plug-in  hybrid  electric,  and  all-electric  vehicles.  Lithium  is  an  ideal  material  for  batteries:  it 
is  the  lightest  metal  in  addition  to  having  a  high  electric  potential.  Charging  a  Li-ion  battery 
usually  means  moving  lithium  cations  from  cathode  into  anode.  There  is  an  increasing 
interest  in  developing  rechargeable  lithium  batteries  with  higher  energy  capacity  and  longer 
cycle  life  for  applications  in  portable  electronic  devices  and  electric  vehicles.  Currently,  the 
graphite  anode  is  the  most  commercially  used  due  to  its  good  capacity  (372  mA  h  g"1) 
together  with  its  rate  capability  and  long  life.  Silicon  has  recently  become  very  popular  as  a 
potential  anode  material  for  lithium  batteries  because  it  has  a  low  discharge  potential  and 
the  highest  known  theoretical  charge  capacity  (which  could  be  lOx  that  of  ghaphite).1  4 
However,  silicon  anodes  have  limited  applications  because  of  the  large  volume  change 
upon  lithium  cations  insertion  or  extraction.  Silicon  nanowires  have  been  shown  to  be 
promising  as  high-performance  lithium  battery  anodes  because  the  can  accommodate  large 
strains  derived  from  lithium  charging  or  discharging. 
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Figure  1.  Si6  skeleton  charges  for  every  global  minimum  structure  vs.  the  number  of 
lithium  atoms  are  plotted.  As  one  can  see,  the  addition  of  five  lithium  atoms  splits  the  Si^ 
skeleton  into  two  fragments:  a  S4  kernel  and  a  two-silicon  fragment. 

Very  recently,  an  ab  initio  study  on  the  lithiation  of  the  Si 4"  cluster  was  reported.  Results  of 
this  study  revealed  that  the  maximum  formal  charge  transfer  from  alkali  metals  to  Si4 
system  is  four  and  that  the  Si 4"  tetrahedral  kernel  is  a  robust  building  block.5  However,  a 
study  on  the  lithiation  of  the  Si6  shows  that  lithium  cations  are  able  to  split  the  silicon 
skeleton  into  two  fragments:  the  S4  kernel  and  a  two-silicon  fragment  (see  Figure  1).  This 
result  actually  gives  an  explanation  to  the  large  volume  change  of  silicon  anodes  upon 
lithium  cations  insertion.  The  latter  study  suggests  that  even  though  the  tetrahedral  SU4' 
structure  is  resistant,  any  material  designed  upon  it  will  be  broken  up  by  lithium  cations.6 
Song  et  al.  have  recently  reported  that  Si/Ge  double-layered  nanotubes  (Si/Ge  DLNT)  are 
promising  materials  as  anodes  for  Li-ion  batteries.  Compared  to  silicon  nanotubes,  Si/Ge 
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DLNT  improve  both  cyclability  and  rate  capability.  It  is  apparent  from  the  last  results  that 
germanium  atoms  play  a  key  role  in  the  silicon  structure  stabilization  upon  lithium  cations 
insertion  and  extraction.  Therefore,  we  have  proposed  to  investigate  both  qualitatively  and 
quantitative  germanium  effects  on  silicon  skeletons. 


Computational  details 

In  order  to  find  out  if  silicon  nanostructured  materials  resist  lithium  cations  insertion  and 
extraction  under  the  presence  of  germanium  atoms,  we  have  been  studying  theoretically  the 
lithiation  process  of  the  SixGey'  (x=4-10,  y=l-10)  clusters.  At  the  moment,  we  have 

explored  the  potential  energy  surfaces  of  the  SixGevLiz'  (x=6,  y=l-2,  z=0-10)  clusters 
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throughout  global  minimum  search  techniques,  i.e.,  Kick  Coalescence  method. 

In  order  to  find  the  global  minimum  structure  of  any  cluster,  a  comprehensive  scanning  on 
the  corresponding  potential  energy  surface  has  been  performed  using  the  Kick  Coalescence 


method.  Kick  Coalescence  (KC)  method,  is  particularly  powerful  given  that  it  combines  the 
accuracy  of  ab  initio  and  DFT  approximations  with  a  very  fast  scanning  of  any  potential 
energy  surface;  in  addition,  it  is  unbiased.  In  this  method,  a  very  large  population  of 
structures  that  are  generated  randomly  is  subjected  to  a  coalescence  procedure  with  all 
atoms  being  pushed  gradually  to  the  molecular  center  of  mass  to  avoid  fragmented 
structures.  Subsequently,  these  structures  are  optimized  to  the  nearest  local  minimum  using 
an  electronic  structure  package.  Currently,  the  KC  method  is  coupled  to  Gaussian  09;  but  it 
can  be  easily  coupled  to  any  software  that  calculates  electronic  structure  of  atoms  and 
molecules.  At  this  stage,  a  low  level  of  theory  is  required  in  order  to  get  a  fast  scanning.  In 
previous  works  on  silicon  clusters,  reliable  results  have  been  got  using  the  combination  of 
the  hybrid  functional  B3LYP  with  the  3-2 1G  basis  set.  This  optimization  process  was 
followed  by  a  reoptimization  and  frequency  calculation  at  a  higher  level  of  theory; 
B3LYP/def2-TZVPP  has  been  used  in  the  present  research. 


Results 

At  the  moment,  we  have  explored  the  potential  energy  surface  of  Si(,Gc",  Si6GeLi5', 
Si6GeLi7',  SieGcLi  m",  Sir,Gc2’,  SieGciLA,  and  SieGcmLi  m"  clusters.  Low-lying  isomers 
structures  are  depicted  in  Figure  2.  As  is  can  be  seen,  Si6  skeleton  is  distorted  in  order  to 
accommodate  the  Ge  atom  in  S^Ge'.  The  global  minimum  structure,  1A,  corresponds  to  a 
planar  five-member  ring  shaped  by  four  Si  atoms  and  one  Ge  atom,  with  other  two  Ge 
atoms  located  above  and  below  the  ring.  The  second  most  stable  isomer,  IB,  is  14.0 
kcal-mol'1  less  stable  than  1A. 
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Figure  2.  Low-lying  isomers  of  Sir,Ge",  Si6GeLi5',  SieGcLi?",  Si^GcLi  |0",  SieGea",  and 
Si6Ge2Lis‘  clusters  calculated  at  the  B3LYP/def2-TZVPP  level.  Si  atoms  are  represented  by 
dark-pink  spheres,  Ge  by  white  spheres,  and  Li  by  light-pink  spheres. 


Given  that  Si6  is  fragmented  when  five  Li  atoms  are  added,  we  decided  to  explore  the 
insertion  of  five  Li  atoms  into  the  S^Ge"  cluster.  The  global  minimum  for  Si6GeLi5' 
corresponds  to  structure  2A,  which  consists  of  a  planar  five-member  Si  ring  surrounded  by 
the  remaining  Si  atoms,  one  Ge  atom  and  five  Li  atoms.  It  is  noticeable  that  Si6  skeleton  is 
not  longer  fragmented  by  the  five  Li  atoms.  However,  there  are  a  number  of  isomers  lying 
in  the  range  of  0  -  1  kcal-mol"1,  where  the  Si6  skeleton  is  actually  fragmented.  These  results 
suggest  that  even  though  one  Ge  atom  helps  to  prevent  Si  fragmentation,  it  is  not  enough. 
In  order  to  corroborate  our  conclusion,  we  decided  to  explore  Sif,GcLi7",  and  SkGcLi  K/.  As 
it  can  be  seen  in  Figure  2,  silicon-germanium  skeleton  is  fragmented  by  Li  atoms  in  both 
cases. 

Subsequently,  we  decided  to  explore  the  effect  of  the  addition  of  two  Ge  atoms  to  S ie 
cluster.  Only  two  low-lying  isomers  were  found:  5A  and  5B.  Isomer  5A  resembles  chair 
conformation  of  cyclohexane  with  the  Ge  atoms  located  above  and  below  the  Si  ring.  The 
second  most  stable  isomer,  5B.  can  be  seen  as  two  Si  tetrahedrons  bridged  by  the  Ge  atoms. 
The  relative  energy  between  both  isomers  is  only  0.6  kcal-mol'1.  When  five  Li  atoms  are 
added  to  the  Si6Ge2_  cluster,  the  resulting  minimum  structure  corresponds  to  6A.  This 
structure  can  be  seen  as  a  silicon  aggregate  surrounded  by  Ge  and  Li  atoms.  However,  the 
second  most  stable  isomer,  6B,  is  a  silicon-fragmented  structure.  Nonetheless,  it  is 
noticeable,  that  the  difference  in  energy  is  2.0  kcal-mol"1  now.  This  result  suggests  that 
incrementing  the  Ge  atoms  into  the  silicon  cluster  actually  prevents  its  fragmentation. 
Because  of  that,  we  decided  to  study  the  biggest  system  proposed,  SiioGeioLiio". 

After  exploring  the  potential  energy  surface  of  SiioGeioLiio",  the  structure  depicted  in 
Figure  3  was  found  to  be  the  global  minimum.  This  structure  accommodates  two  distorted 
silicon  tetrahedra  connected  by  a  silicon  atom,  surrounded  by  other  three  silicon  atoms  as 
well  as  the  germanium  fragment  divided  into  two  skeletons  that  also  corresponds  to  two 
distorted  tetrahedra.  The  lithium  atoms  are  accommodated  around  the  Si-Ge  skeleton. 
There  is  no  fragmentation  due  to  lithium  atoms  in  this  structure.  The  last  result  confirms 
that  actually  germanium  atoms  prevent  fragmentation  of  the  silicon  structure,  as  well  as  the 
tetrahedral  skeleton  of  silicon  atoms  can  be  used  as  building-block  fragment  that  stabilizes 
these  materials. 
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Figure  3.  Global  minimum  structure  of  SiioGeioLiio’  (7A). 


The  future  work  consists  of  studying  materials  based  on  the  tetrahedral  skeleton  of  silicon 
atoms,  surrounded  by  tetrahedral  fragments  of  germanium  atoms.  Such  a  material  should 
resist  fragmentation  and  be  adequate  to  be  used  as  anode  in  Li-ion  batteries. 

Finally,  the  xyz  Cartesian  coordinates  are  described  below  in  order  to  provide  better  detail 
of  the  studies  in  this  work.  These  Cartesian  coordinates  can  be  displayed  in  programs  such 
as  Mercury,  ArgusLab,  etc,  or  some  other  program  that  supports  this  type  of  coordinates. 
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